Phase instability and coarsening in two dimensions 
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Instabilities and pattern formation is the rule in nonequilibrium systems. Selection of a persistent 
lengthscale, or coarsening (increase of the lengthscale with time) are the two major alternatives. 
When and under which conditions one dynamics prevails over the other is a longstanding problem, 
particularly beyond one dimension. It is shown that the challenge can be defied in two dimen- 
sions, using the concept of phase diffusion equation. We find that coarsening is related to the 
A — dependence of a suitable phase diffusion coefficient, Dn(A), depending on lattice symmetry and 
conservation laws. These results are exemplified analytically on prototypical nonlinear equations. 

PACS numbers: 05.70.Ln, 05.45.-a, 64.60.Ht, 68.43. Jk 
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Introduction. — Pattern formation (or morphogenesis) 
is abundant in nature, both in inanimate and living sys- 
tems. Examples [l[ are encountered in many branches of 
science: physics (e.g. sand ripples), chemistry (chemical 
spots, reminiscent of those on animal skins, e.g. jaguar), 
biology (asters of macromolecules during cell division), 
and so on. Patterns arise often due to the loss of stabil- 
ity of an initially structureless state. The diversity and 
richness of morphogenesis is concomitant to the nonlinear 
and nonequilibrium nature of these systems. 

A classification of dynamics and patterns that prevail 
for a given nonlinear system is a challenging task. For ex- 
ample, given a system described by non linear equations, 
it is not obvious to state a priori whether one would assist 
to the selection of a pattern with a specific lengthscale 
during time, or rather coarsening (increase of lengthscale 
with time) would prevail. 

Coarsening is also the primary scenario in phase sep- 
aration processes [2j and it attracts a continuous inter- 
est 0, i, H, H, SIS H GU motivated by the demand 
for a deeper understanding. This goal may be attained 
following two main directions: either investigating rig- 
orous solutions to specific models 0, Q, or developing 
general approaches which are valid for wide classes of 
models/equations. Regarding the second strategy, the 
only example we are aware of is due to A. Bray and col- 
laborators 0, [ll| , who analyzed some models which are 
derivable from a potential: comparing the global rate of 
energy change with the energy dissipation, the temporal 
growth law for the typical lengthscale, X(t), could be de- 
rived. The present study can be seen as complementary, 
in two respects. First, we do not impose to the equa- 
tions to be derivable from a potential. Second, we do 
not assume that coarsening does occur, we are able to 
determine if and when it occurs, instead. Note that our 
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study does not address the other kind of coarsening which 
consists in an increase of a domain with a certain orien- 
tation (or pattern) , and / or lengthscale, at the expense of 
another domain [T3 |. 

A few years ago we established a coarsening criterion 
in one dimension (ID) Q, for certain classes of nonlin- 
ear equations, generalizations of the celebrated models A 
and B of the dynamics [l3| . The criterion demands to 
analyse the stationary periodic states, which are found 
to solve Newton's equations where the spatial variable 
plays the role of time and the fictitious particle oscillates 
in a potential well of arbitrary shape (different equations 
give rise to different potentials). The oscillation period 
of the particle corresponds, in the mechanical analogy, 
to the wavelength A of the periodic steady state and it 
is a function of the amplitude A, A = X(A). The crite- 
rion simply states that coarsening occurs if and only if 
dX/dA > 0. 

Our analysis anticipated the existence of an "inter- 
rupted coarsening" scenario, which has recently been in- 
voked fl4| in the context of wind driven sand dunes. It 
is clear that these physical phenomena (as well as, e.g ., 
mound coarsening at the nanoscale in crystal growth [151 ]) 
demand for the construction of a full two-dimensional ap- 
proach, which is the basic goal of the present paper. 

However, in our former ID criterium, the spatial vari- 
able plays the role of time (see above). Therefore, it is 
obvious that this concept is limited to one spatial dimen- 
sional systems and extensions of the criterion to higher 
dimensions seemed to present a significant challenge. It 
is shown here that the challenge to find general crite- 
ria for understanding if a two dimensional (2D) model 
shows coarsening or not can be defied, following the con- 
cept of phase diffusion equation [l6j|. Before discussing 
the details, we give a simple overview of the method and 
results. 

The method. — A perfectly periodic steady state is de- 
fined through a pair of wave wectors qi , q2 and it is a 
function of the phases 4>i,2 = qi,2 ' x - Perturbations 
make qi ; 2 acquire slow dependencies on time and space, 
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= q;(X, T). The dynamical response to perturbations 
is described by phase equations, which have, at the lowest 
order, the form of linear diffusion equations: 
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Setting 0i 2 = <f>i 2 exp(ajT + iK ■ X), standard Fourier 
analysis allows to determine the stable/unstable charac- 
ter of phase dynamics. Coarsening is related to instabil- 
ity, which is signaled by a positive oj. The requirement 
lj > for some K implies a condition on the diffusion 
coefficients Djy , which are functions of steady state prop- 
erties only. The central result is that the condition uj > 
finally writes as dA/dX > 0, where A is some function of 
the stationary periodic solution having wavelength A. 

In one dimension, A is shown to correspond to the am- 
plitude of the stationary solution Q. In two dimensions, 
A takes different expressions, depending on the underly- 
ing symmetry of the periodic stationary solutions, and it 
dees do not seem to have a simple physical interpreta- 
tion. Nevertheless, the condition dA/dX > docs have 
a simple reading: coarsening occurs if and only if A is 
an increasing function of A. That is to say the coarsen- 
ing criterion can still be established only upon inspecting 
steady-state solutions. Even more importantly, the coars- 
ening law, X(t), which describes the time dependence of 
the typical size of the pattern, is obtained through a re- 
lation of the form D(X)t ~ A 2 (or equivalcntly to ~ 
where D is a typical diffusion constant. In the limit of 
large time, a power law behavior, A ~ t n , is expected. 

We shall exemplify the method on classical model 
equations. For instance, we show that n — \ f° r the 
non-conserved real Ginzburg-Landau equation (model 
A), and n — g for the conserved Cahn-Hilliard equation 
(model B). As will be recognized, the interesting mes- 
sage is that the methodology does not evoke whether the 
model equation are variational or not. Thus the study 
should work with any other equation, be it of potential 
(i.e. variational) nature or not. 

Generalized Ginzburg-Landau type models. — We now 
apply and discuss our method for a nonconserved class 
of equations, 



d t u = A(u) + V 2 u, 



(2) 



which reduces, for A(u) = u — u 3 , to the famous real 
Ginzburg-Landau (or Allen-Cahn) equation. This equa- 
tion is known to exhibit perpetual coarsening, with n = 
\ 0] . Here we only require A(u) « u for small u, other- 
wise A(u) can be any function of u. 

A steady periodic state Uo(x) depends on the fast spa- 
tial variables x = (x\,X2), or equivalcntly 0x,2 = Qi,2 ■ x 
which are the fast (constant) phases. Phase modes are 
studied by perturbing the perfectly periodic steady state. 
The perturbations of interest depend on slow spatial and 
temporal variables. The slow character of relevant per- 
turbations reflects the property of the Goldstone mode (if 
m (x) is a periodic solution, so is u (x + Ro)), which has 
infinite relaxation time and is a quasi-dangerous mode, 
making slowly varying perturbations to be most persis- 
tent. Let us encode the slow modulation by a dimension- 
less parameter e <C 1. We define slow spatial and time 
scales as 



Xi = exi, X 2 — ex 2 , T 



e 2 t, 



(3) 



where the factor e 2 is an ansatz dictated by the fact that 
phase rearrangement occurs via diffusion. We introduce 
the slow phases 1/11,2 = e0i,2i so that = V x <fo = VxV'i) 
where X = {X^ X-z) = ex. 

In general terms, at order e we can write Q 



\dipi dip 2 
~ e ~dT 01 ~dT 02 
V x = qi5«k + q 2 <9 lfc + eV x 



(4) 
(5) 



The above expressions for the derivative, along with 
the standard expansion u = uo + eux, are substituted into 
Eq. ([2]) . The zeroth order problem yields the differential 
equation obeyed by the two-dimensional profile of steady 
states, 



Afluo] = A(u ) + V 2 u = 



(6) 



where the subscript in the V 2 operator is taken to mean 
that it is evaluated for e = in Eq. ((5|). 
The next order in e has the form na, [l3l 



C[ui] = g{ua,d T il>i,dx i x j il)k), 



(7) 



where C = Co = (A'(uq) + V 2 ) is the linear operator 
obtained as Frechet derivative of J\f, and 



9 = + ~§F^ 2 ~ t 2 ^ 1 ^ 1 + q2 ^) ' V x wo + (V x • qi)9^Mo + (Vx ■ 02)^2^0] • (8) 



The expression for g can be rewritten as where Vi = dfcUo, qki is the i— th component of q/c and 

g = (d T ipi)vi + (drip2)v2 - 2q ki (ipi) ij - (ipi)uVi (9) 

oqij 
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Since £q = Cq for Eq. ([2|), wc easily get u>, = Vi. 
leads to the sought-after phase diffusion equations, 



=d 2 4 H /dx l dx 3 - 

In order to avoid secular terms in Eq. (J7J), if the homo- 
geneous equation £ [w] = has a nonvanishing solution, 
w must be orthogonal to g [IH, (w,g) — (the precise 
definition of the inner product is given below). 

There exit two nonvanishing solutions w\ , w 2 in two 
dimensions, resulting thus in two solvability conditions. 

This 

(10) 
(11) 

(12) 
(13) 



{v 2 )dTtpi + (viV 2 )dT4'2 



E 
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i J 



where (v 2 ) = (vf) = (v 2 ) and 
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k ql] 

^2lki{v2^^) + Sij(v 2 Vi). 
k 



There are five Bravais lattices in 2D (oblique, rect- 
angular, centered rectangular, hexagonal, and square). 
We have performed explicitly the calculation for squares, 
hexagons and rectangles. Here, wc shall focus on the 
first two symmetries. For square symmetry, qi = g(l,0), 
q2 = <z(0, 1) and (viv 2 ) = 0. Finally, we get 



dT 

dip 2 
dT 



d 2 j; 2 

dx 2 
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dX x X 2 
dX t X 2 



(14) 



with Dn = D\\ = d q (q(v 2 ))/(v 2 ) and L>i 2 = D\l = 
D\\ — l. A linear stability analysis is performed by setting 
i/>i,2 = 2 exp(wT + iK ■ X) in Eqs. (14]), which yields 



■w(l + D n )K 2 + D n K 4 = 



(15) 



whose roots are u>\ = —K 2 and lo 2 = —DnK 2 . 

While oj\ < 0, the sign of c<j nc = uj 2 depends on the sign 
of d q A = d g (q{v 2 )). Therefore, oj 2 is positive, implying 
phase instability, if and only if ^4. = q(v 2 ) is an increasing 
function of the wavelength A. 

For hexagonal symmetry, qi i2 = q{\,±^) and ^1,2 
must be linearly combined in order to get Eqs. (|14l) . with 
a different 



D n = D\T 



d q (q(v 2 ) 2 ) 
(v 2 ) 2 



(16) 



albeit the relation D\ 2 X = f ff — 1 still holds. 

Cahn- Hilliar d type models. — Let us now consider a con- 
served equation, a generalized form of the well-known 
Cahn-Hilliard equation, 



d t u = -V 2 (A(u) + V 2 



(17) 





square 


other symmetries 


standard GL 


1/3 


1/2 


modified GL 


1/2 


1/2 


standard CH 


1/4 


1/3 



TABLE I: Summary of coarsening exponents for different 
models and base symmetries. 



We only provide the result. Wc formally obtain 
Eq. {?]), but now C = — Vg£o and g = {dr^ijvi + 
(dTip 2 )v 2 + VqV 2 u where V 2 u is a shorthand nota- 
tion for the expression in square parentheses in Eq. |8|). 
C is not self-adjoint, because & = — £qVq and its ker- 



nel 



so that 



is such that C^wi = = — £oVqc 

- — d"o 
w i — v * — d<t>i ■ 

The orthogonality conditions (uiig) — are a bit more 
challenging to analyse, and will be discussed elsewhere. 
Suffice it here to say that the relevant eigenvalue reads 



(18) 



where lo c and w NC are the eigenvalues associated with 
the unstable mode for the Conserved and Non Conserved 
versions of a given model, and I c = (v 2 )/(wv) is the 
factor taking into account the conservation law. 

The coarsening exponent. — The coarsening law is de- 
termined on the basis of the relation \uj(q)\ ~ 1/t, where 
u) = w NC = —D\\q 2 or oj = lo c = I c lo nc . A major role is 
played by the quantity 



1 



r2x p2n 

Jo #1 /o ^ 



duo 



1 



(2tt) 2 



(19) 



Consider first the Ginzburg-Landau (nonconserved) 
equation. In this case, uo is a constant except along do- 
main walls, so (v 2 ) scales linearly with A. More precisely, 
(v ) = a/q + b. Reporting this expression into Dn yields 
n = i for square symmetry and n = i for hexagonal 
symmetry. Preliminary analysis for rectangular symme- 
try leads to n — i as well. In general had (v 2 ) scaled 
as l/q a (a > 1) wc would have obtained n = | for all 
symmetries. Therefore, GL scaling ((v 2 ) ~ 1/q) obtained 
for the square symmetry seems to be singular: all other 
cases provide n = \. 

The computation of the coarsening exponent for con- 
served models requires estimation of I c = (v 2 )/(wv). 
Since V 2 ui = v, we infer (after integration by parts) 
I c ~ {v 2 )/( u o) ~ 1- Therefore, for the standard Cahn- 
Hilliard model, we get n = ^ for square symmetry 
and n = \ for other symmetries. Table U summarizes 



the results. Note that standard GL and CH refer to 
A(u) = 



u , so that {v 
„2\ _ 



q . In contrast, modified 
GL equation means (v 2 ) ~ q" a with a > 1. 

Discussion. — Wc have found that the 2D coarsening 
exponent for the standard GL and CH equations is the 
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expected one (n = | and n = |, respectively) for all 
symmetries but the square one, which exhibits a slower 
coarsening . The peculiar behaviour of square symmetry 
is at present not understood. It must be noted, however, 
that if (v 2 ) ~ q~ a (a > 1, as occurs with the modi- 
fied GL) then n = \ for any symmetry. It seems thus 
that the peculiar behavior of square symmetry is quite 
specific to the considered equation, rather than general. 
Nonetheless, a deeper understanding of this fact merits 
higher attention in the future. 

It is worthstressing once again that our approach al- 
lows to reduce the study of dynamical properties (coars- 
ening) to the behavior of quantities (flu, Dn) depending 
on steady states properties only. Our analysis has been, 
for ease of presentation, exemplified on the two classical 
equations, namely the GL and CH ones, but the method 
can be applied to other equations. In particular, our ap- 
proach (as it is evident in Ref. [l7| ) does not require that 
the equation is derivable from a potential. 

While for the GL and CH equations we could extract 
analytically the coarsening exponent, it may prove nec- 
essary that for other equations there would be a need for 
numerical solutions of the steady-states problem in order 
to evaluate the diffusion constants. This is a quite simple 
task numerically, even at higher dimensions (where full 
time-dependent studies are difficult, or even unfeasible if 
the asymptotic regime is to be ascertained). 

The amplitude branch is stable, in contrast to the 
phase one when coarsening is to take place. There are 
five different branches corresponding to different symme- 
tries (five Bravais lattices) . The coarsening of the square 
branch is slower (only for the classical GL and CH equa- 
tions) than that of the other symmetries. It is natural to 
expect the fastest growing structure to prevail. Symme- 
tries other than the square one offer, for these classical 
equations, a faster channel towards coarsening. 

A particularly important question concerns relevance 



of slower coarsening peculiar to square symmetry. A 
possibility might be offered [Til] by the Rayleigh-Benard 
convection, which may develop cither rolls or hexagons, 
with different local symmetries in different spatial re- 
gions. These regions coarsen in time and the coarsening 
rate might be different for domains of rolls and hexagons. 

We have investigated steady-state 'periodic solutions, 
while during coarsening, even in ID, no long range or- 
der is observed. The questions thus arises about impor- 
tance of periodic solutions. Locally in space the ampli- 
tude is adiabatically slaved to the phase, so that (for a 
given wave vector) it reaches the quasi-steady solution. If 
steady-state solutions were stable (no coarsening), then 
the system would generically choose the periodic solu- 
tion (possibly with defects), as is known for many pat- 
tern forming systems. Thus, having shown here that the 
periodic solutions are unstable with respect to phase fluc- 
tuations, we expect that the pattern should coarsen. One 
can not exclude, however, the existence of non-periodic 
(like disordered) stable solutions for nonlinear extended 
systems with an average wavelength which does not grow 
in time. We are not aware of any such scenario, however. 

Finally, it is noteworthy that the phase diffusion equa- 
tions share the same structures as those for the dis- 
placement field of 2D cry stals [20I ] (see also the Sup- 
plementary Materials |21(). Conservation law imposes 
dt<j>i = —dkJik (this is the analogue of the dynamical 
equation in continuum media) where Jik is the phase 
current. Jik is proportional to the gradient of the phase, 
di<fi m , and the proportionality coefficient is a rank four 
tensor, Jik = Aku m di4> m (this is the analogue of Hooke's 
law). A should be invariant under the symmetry group 
of the considered crystal. The phase equations has thus 
the same structure as the elasticity problem of crystals. 
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Supplementary materials 

Derivation of the phase equation from symmetry 
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Like in crystal elasticity, the phase equation has the form 



dtipi = -d k Jik 

where Jik is the phase current which reads 

J lk = -A lk i m (diip m + d m ipi). 

This is the analogue of Hooke's law. The minus sign in front A expresses the fact that the current is opposite to the 
gradient. Here we have written the current in a symmetrized form, since like in elasticity, the deformation tensor 
(which is a measure of distances in elasticity) is symmetric (the antisymmetric part corresponds to global rotation of 
the pattern, and is thus unimportant). 
Thus the phase equation becomes 

dtipi = -Aikimdkiditpm + 9 m ipi) , (20) 

where A, which is the analogue of the matrix defining the Lame coefficients, should be invariant under the symmetry 
group of the considered crystal. Let us denote differentiation with respect to the variables X\ and X2 simply by d\ 
and 82, respectively. In a cubic crystal, as well as in a square lattice, only three components of A are non zero (see 
Rcf.[l]), 

(21) 
(22) 
(23) 

It follows that 

dtipi = 2X 1 dx 1 ip 1 + (2A 2 + X 3 )d 12 ip2 + Xsd 2 %ipi , (24) 

which has the form of our first equation (14). The same reasoning leads to the second equation (14). Regarding the 
first equation, term 822^1 is absent from symmetry, while in the second equation, dntpi is absent. 
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